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^ Abstract 

We consider a nonlinear 4th-order degenerate parabolic par- 
tial differential equation that arises in modelling the dynamics 
of an incompressible thin liquid film on the outer surface of 

1 a rotating horizontal cylinder in the presence of gravity. The 

parameters involved determine a rich variety of qualitatively 

^ different flows. Depending on the initial data and the parame- 

^ ter values, we prove the existence of nonnegative periodic weak 

^ solutions. In addition, we prove that these solutions and their 

gradients cannot grow any faster than linearly in time; there 
^ cannot be a finite-time blow-up. Finally, we present numerical 

^ simulations of solutions. 

as 
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1 Introduction 

We consider the dynamics of a viscous incompressible fluid on the 
outer surface of a horizontal circular cylinder that is rotating around 
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Figure 1: Liquid film on the outer surface of a rotating horizontal 
cylinder in the presence of gravity. 

its axis in the presence of gravity, see Figure [T| If the cylinder is fully 
coated there is only one free boundary: where the hquid meets the 
surrounding air. Otherwise, there is also a free boundary (or contact 
line) where the air and liquid meet the cyhnder's surface. 

The motion of the liquid film is governed by four physical effects: 
viscosity, gravity, surface tension, and centrifugal forces. These are 
reflected in the parameters: R — the radius of the cylinder, uj — its rate 
of rotation (assumed constant), g — the acceleration due to gravity, u 
— the kinematic viscosity, p — the fluid's density, and a — the surface 
tension. 

These parameters yield three independent dimensionless numbers: 
the Reynolds number Re = {R^uj)/^, the Galfleo number Ga = (7/ {Ruo"^) 
and the Weber number We = {pR^uj'^)/a. 

We introduce the parameter e = h/R^ where h is the average thick- 
ness of the hquid. The following quantities are assumed to have finite, 
nonzero hmits as e ^ 0. [2S1 [211 [21 121: 

K = Ree^, X = =^e^, and /i = GaRee^. (1.1) 

This corresponds to a low rotation rate, for example. 

One can model the fiow using the full three-dimensional Navier- 
Stokes equations with free boundaries: for m(x, ?/, 2;, t) in the region 
X e [— tTjTt), y G and z G {0, h{x,y,t)) where x is the angular 
variable, y is the axial variable, and h{x,y,t) is the thickness of the 
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fluid above the point (x, y) on the surface of the cyhnder at time t. 
This has been done by Pukhnachov |28) in which he considered the 
physical regime for which the ratio of the free-fall acceleration and the 
centripetal acceleration is small. There, he proved the existence and 
uniqueness of fully-coating steady states (no contact hne is present). 
We know of no results for the affiliated initial value problem. 

In this physical regime, if one also makes a longwave approxima- 
tion (the thickness of the coating fluid is smaUer than the radius of the 
cyhnder) and if one further assumes that the rotation rate is low (or the 
viscosity is large) then the three-dimensional Navier-Stokes equations 
with free boundary can be approximated by a fourth-order degener- 
ate partial differential equation (PDE) for the fllm thickness h{x,y,t). 
This is done by averaging the fluid flow in the direction normal to the 
cyhnder |28l [2^. If one further assumes that the flow is independent 
of the axial variable, y, then this results in a PDE in one dimension 
for h{x, t). 

In his pioneering 1977 article about syrup rings on a rotating roller, 
Moffatt neglected the effect of surface tension (i.e. We~^ = = x), 
assumed the flow was uniform in the axial variable, and derived [21] 
the following model for the thin film thickness: 

ht + {h- !lh^ cos{x))^ = 0, (1.2) 

where fi is given in ( 1.1[ ) and 



X G [— 7r,7r], t > 0, h is 27r-periodic in x. 

Pukhnachov's 1977 article [28j gives the first model that takes into 
account surface tension: 

ht + ih- Ih^ cos(x)), + f {h' + K,,))^ = (1.3) 



where /i and x ctre given in (1.1) and 

X e [— 7r,7r], t > 0, h is 27r-periodic in x. 
This model assumes a no-slip boundary condition at the liquid/sohd in- 



terface. For a solution to (1.2) or (1.3) to be physically relevant, either 
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h is strictly positive (the cylinder is fully coated) or h is nonnegative 
(the cylinder is wet in some region and dry in others). 

Surprisingly little is understood about the initial value problem for 



(1.3). Bernis and Friedman were the first to prove the existence of 
nonnegative weak solutions for nonnegative initial data for the related 
fourth-order nonlinear degenerate parabolic PDE 

ht + U{h)h...)x = ^, (1.4) 

where /(/i) = |/ir/o(/i), /o(/i) > 0, n^l. 

Unlike for second-order parabolic equations, there is no comparison 



principle for equation (1.4). Nonnegative initial data does not auto- 
matically yield a nonnegative solution; indeed it may not even be true 
for general fourth-order PDE (e.g. consider ht = —h^xxx)- The de- 



generacy f{h) in equation (1.4) is key in ensuring that nonnegative 
solutions exist. 



Lower-order terms can be added to equation (1.4) to model addi- 
tional physical effects. For example, 

ht + {f{h) hxxx)x - {g{h)hx)x = (1.5) 



where g{h) > for h ^ 0. Equation (1.5) can model a thin liquid film 
on a horizontal surface with gravity acting towards the surface. If this 
surface is not horizontal then the dynamics can be modelled by 

ht + {h''{a-bK + Kxx))x = 0, a>0, 6>0 (1.6) 

The constant a in the first-order term vanishes as the surface becomes 
more and more horizontal. If the thin film of liquid is on a horizontal 
surface with gravity acting away from the surface then the thin film 
dynamics can be modelled by 

ht + {f{h) hxxx)x + {g{h)hx)x = 0. (1.7) 
For a thorough review of the modelling of thin liquid films, see [1 



In equations (1.5) and (1.6) the second-order term is stabilizing: 



if one hnearizes the equation about a constant, positive steady state 
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then the presence of the second-order term increases how quickly per- 



turbations decay in time. In equation (1.7), the second-order term is 



destabilizing: the linearized equation can have some long-wavelength 
perturbations that grow in time. For this reason, we refer to equation 



(1.7) as "long- wave unstable". The long-wave stable equations (1.5) 



and (1.6) have similar dynamics as equation (1.4) however the long- 



wave unstable equation (1.7) can have nontrivial exact solutions and 



can have finite-time blow-up {h{x*,t) | oo as t 1 1* < oo). 

In all cases, the fourth-order term makes it harder to prove desir- 
able properties such as: the short-time (or long-time) existence of 
nonnegative solutions given nonnegative initial data, compactly sup- 
ported initial data yielding compactly supported solutions (finite speed 
of propagation), and uniqueness. Indeed, there are counterexamples 
to uniqueness of weak solutions [3|. Results about existence and long- 
time behavior for solutions of (1.5) can be found in |6]; analogous 



results for (1.6) are in [TH]. See [HIS] for results about existence, finite 



speed of propagation, and finite-time blow-up for equation (1.7). 



In this paper we study the existence of weak solutions of the thin 
film equation 

ht + {\h\^{ao h^^^ + aih^ + a2w'{x)))^ + Og/i^ = (1.8) 
where ai, 02, 03 are arbitrary constants, constant Oq > 0, and w{x) 



is periodic. Equation (1.3) is a special case of (1.8). The sign of ai 
determines whether equation (1.8) is long-wave unstable. Also, the 



coefficient of the convection term a2{w' {x)\h\ )x can depend on space 
and will change sign if a2w'{x) ^ 0. The cubic nonhnearity \h\^ in 



equation (1.8) arises naturally in models of thin liquid films with no- 
shp boundary conditions at the liquid/sohd interface. Our methods 
generahze naturally to f{h) = we refer the reader to |2l El E] for 
the types of results expected. 

Given nonnegative initial data that satisfies some reasonable condi- 
tions, we prove long-time existence of nonnegative periodic generalized 



weak solutions to the initial value problem for equation (1.8). We start 



by using energy methods to prove short-time existence of a weak solu- 
tion and find an explicit lower bound on the time of existence. A gen- 
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eralization and sharpening of the method used in |H] allows us to prove 
that the norm of the constructed solution can grow at most linearly 
in time, precluding the possibility of a finite-time blow-up. This 
control, combined with the explicit lower bound on the (short) time of 
existence, allows us to continue the weak solution in time, extending 
the short-time result to a long-time result. 



If 02 = or 03 = in equation (1.8) then solutions will be uniformly 



bounded for all time. If 02 7^ and 03 7^ 0, it is natural to ask if the 
nonhnear advection term could cause finite-time blow-up {h{x*,t) | 00 
as t t t*. Such finite-time blow-up is impossible by the linear-in-time 
bound on but we have not ruled out that a solution might grow in 
an unbounded manner as time goes to infinity. 

In [m [15], the authors consider the multidimensional analogue of 



(1.4) 



+ V- (|/irVA/i) = 0, (1.9) 

for h{x,t) where x G f2 C with N = 2,3. Depending on the sign 
of A', if g = then equation 

ht + V- {f{h)VAh + VA{h)) = g{t, x, h, Vh) (1.10) 



on is the multidimensional analogue of equation ( 1.5[ ) or (1.7). In 



[T6] . the authors consider the long- wave stable case with g = and 
power-law coefficients, f{h) = \h\^ and A'{h) = —\h\'^. In [I9|, the 
author considers the Neumann problem for both the long-wave stable 
and unstable cases with the assumption that f{h) > has power-law- 
like behavior near h = 0, that |A'(/i)| is dominated by f{h) (specifically 
|y4'(/i)| ^ dof{h) for some do), and that the source/sink term g{t,x, h) 
grows no faster than hnearly in h. In [33l|3ll|36], the authors consider 



the Neumann problem for the long- wave stable case of (1.10) with 
power-law coefficients and a larger class of source terms: g{t, x, h) ~ 
\h\^~^h with. \ > 0. In |3n[35]. the same authors consider the long- wave 
stable equation with power-law coefficients but with g{h) = a ■ Vh{h) 
where h{z) ~ and a G M^: g models advective effects. They consider 
the problem both on and on a bounded domain f2. 



All of these works on (1.9) and (1.10) construct nonnegative weak 



solutions from nonnegative initial data and address qualitative ques- 
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tions such as dependence on exponents n and m and A, on dimension 
N, speed of propagation of the support and of perturbations, exact 
asymptotics of the motion of the support, and positivity properties. 
We note that the works [331 EH [Ml EH ES] also construct "strong" 
solutions. 

Finally, we refer readers to the technical report [13] which presents 
the results of this article, and some additional results, along with more 
extensive discussion, calculations, and simulations. 



2 Steady state solutions 



Smooth steady state solutions, h{x,t) = h{x), of (1.3) satisfy 



h - f /i^ cos(a;) + f [Y? ih^ + h.,^^)) = q (2.1) 

where g is a constant of integration that corresponds to the dimension- 
less mass flux. In the zero surface tension case (x = 0), steady states 
satisfy 

h - cos(x) = q. (2.2) 

Such steady states were first studied by Johnson [20] and Moffatt [24j. 
Johnson proved that there are positive, unique, smooth steady states 
if and only if the flux is not too large: < g < 2/(3y^). These steady- 
states are neutrally stable [26]. Smooth, positive steady states in the 
presence of surface tension have been studied by a number of authors. 
One striking computational result [T] is that for certain values of x and 
/i there can be non-uniqueness. 

These non-unique steady states were numerically discovered via an 
elegant combination of asymptotics and a two-parameter (mass and 
flux) continuation method [H Figure 14]. To start the continuation 
method, earlier work [2] on the regime in which viscous forces dominate 
gravity was used. There, asymptotics show that for smaU fluxes the 
steady state is close to g + l/3g^ cos(x) + O(g^), providing a good flrst 
guess for the iteration used to flnd the steady state. The bifurcation 
diagram shown in Figure 14 of [1] also suggests that the Moffatt model 
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(1.2) can be considered as the limit of the Pukhnachov model (1.3) as 
surface tension goes to zero (x — * 0). 

Pukhnachov proved [30] a nonexistence result: no positive steady 
states exist if g > 2^/^JJ^ ~ 3.464/^//!. We improve this, proving that 
no such solution exists if g > 2/3 A/27/i ~ 0.943/^/1. 

Proposition 2.1. There does not exist a strictly positive 2n periodic 



solution h{x) of equation (2.1) if q > 2/3 \/2/ij,. 

Proof of Proposition \2. 1\ Following Pukhnachov, we start by rescal- 
ing the flux to 1 by introducing y{x) = h{ x)/q and introducing the 

3 2 I 

parameters 7 = ^ and P = Equation (2.1) transforms to 



l{y"' + y')=f3cos{x)-^ + ^. (2.3) 

The solution y is written as y{x) = Oq + ai cos(a;) + a2sin(a;) + v{x) 
where v{x) _L span{l, cos(x), sin(a:)} and satisfies 

^{v"' + v')=Pcos{x)-^ + ^. (2.4) 



A solution V exists only if the right-hand side of (2.4) is orthogonal to 
span{l, cos(x), sin(x)}. As a result. 



1 \J™_n / / 1 



dx = 



cos(x) dx = 7r/5. (2.5) 



It foUows from (|2^ that vr/? < / ^ (1 + cos(x)) dx < ^ 27r. 



y>l 

This shows that if there is a positive steady state then (3 < 8/27. 
Recalhng the definition of (3, there is no steady state if g > 2/3 ^/2/j2. 

□ 

The proof also holds in the case of zero surface tension x = 7 = 
and so it is natural that the bound 2/3 \/2/n is larger than 2/(3^/1) 
(the bound found by Johnson and Moffatt.) Also, we note that numer- 
ical simulations that suggest nonexistence of a positive steady state if 
q > 0.854 when /i = 1 for a large range of surface tension values 
p. 61]; our bound of 0.943 is not too far off from this. 
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3 Short-time Existence and Regularity of 
Solutions 

We are interested in the existence of nonnegative generalized weak 
solutions to the following initial-boundary value problem: 

{ht + {f{h){aoh^^^ + aih^ + a2w'{x)))^ + a^h^ = in Qt, (3.1) 
g(-a,t) = |$(a,t) for t > 0, 2 = M, (3.2) 
/i(x,0) = /io(a;) ^ 0, (3.3) 

where f{h) = \h\^, h = h{x,t), Q = {-a, a), and Qt = Q x (0,T). 
Note that rather than considering the interval (—a, a) with boundary 



conditions (3.2) one can equally well consider the problem on the circle 
S^; our methods and results would apply here too. Recall that ai, 02, 



and 03 in equation (3.1) are arbitrary constants; Oq is required to be 



positive. The function w in (3.1) is assumed to satisfy: 



w e C^+"'{n) for some < 7 < 1, f^(-a) = ff (a) for z = 0, 2. 

(3.4) 

We consider a generalized weak solution in the following sense [Sill]: 

Definition 3.1. A generalized weak solution of problem (P) is a func- 
tion h satisfying 

h G C'Jf''^\QT) n L-(0,T; H\Q)), (3.5) 
hteL\0,T;{H\Q)y), (3.6) 

h e C^iC?), ^/Jih) (aoh,,, + a^K + a2w') E 1^7), (3.7) 



where T = Qt \ {h = Ut = 0) and h satisfies (3.1) in the following 
sense: 

T 

{ht{-,t),(f)) dt - j j f{h){aoh^^^ + aih^ + a2w'{x))(f)xdxdt 
y 

— 03 / / h(j)x dxdt = (3.8) 
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for all (f> G C^{Qt) with (p{—a, ■) = 0(a, ■); 

h{-,t) h{-,0) = liQ pointwise & strongly in L^(fi) as t 0, (3.9) 
h{-a,t) = h{a,t) Vt e [0,T] and g(-a,t) = g(a,t) (3.10) 
/or i = 1,3 at all points of the lateral boundary where {h ^ 0}. 



Because the second term of (3.8) has an integral over T rather 



than over Qy, the generahzed weak solution is "weaker" than a stan- 



dard weak solution. Also note that the first term of (3.8) uses ht G 
L'^{0,T; (iJ^(n))'); this is different from the definition of weak solution 
first introduced by Bernis and Friedman [5]; there, the first term was 
the integral of kept integrated over Qt- 

We first prove the short-time existence of a generalized weak solution 
and then prove that it can have additional regularity. In Section |4] we 
prove additional control for the norm which then allows us to prove 
long-time existence. 

Theorem 1 (Existence). Let the nonnegative initial data Hq G H^{Q) 
satisfy 

dx < oo, (3-11) 



/ 



ho{x) 



n 

and either 1) ho{—a) = ho{a) = or 2) ho{—a) = ho{a) ^ and 
^^(— a) = 7^(a) holds for i = 1,3. Then for some time Tioc > 
there exists a nonnegative generalized weak solution, h, on Qti„^ in the 



sense of the definition \3.1\ Furthermore, 

heL\0,Tioc;H\n)). (3.12) 

Let 

Eo{T) := ^ J (aohl{x,T) — aih'^{x,T) — 2a2w{x)h{x,T)) dx (3.13) 

then the weak solution satisfies 
^o{Tioc)+ j j h^iaoh^^^ + a^h^ + a2w'fdxdt ^ 8,o{0)+K Ti^c, (3.14) 

{h>0} 
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where K = \a2a3\ \\w'\\ooC < oo. The time of existence, Tioc, is deter- 
mined by ao, ai, a2, w' , \VL\, and Hq. 

We note that the analogue of Theorem 4.2 in [H] also holds: there 
exists a nonnegative weak solution with the integral formulation 

0. ^> * + a„ lj(3h%.k„^. + d.dt (3.15) 

Qt 

— J J {aihx + a2w' + a^h) (p^ dxdt = 0. 

Qt 

Theorem 2 (Regularity). If the initial data from Theorem^ also sat- 
isfies 

J hQ~^{x) dx < oo 
n 

for some —l/2<a<l, a ^ then there exists < T^^J < Tioc such 
that the nonnegative generalized weak solution from Theorem\I\has the 
extra regularity 

h'^ G L\0,Tlf,H''in)) and h"^ e ^^(0, t/^^; 

The solutions from Theorem [2] are often called "strong" solutions in 
the thin film literature. If the initial data satisfies / Hq'^ dx < oo then 
the added regularity from Theorem [2] allows one to prove the existence 
of nonnegative solutions with an integral formulation [7] that is similar 



to that of (3.15) except that the second integral is replaced by the 
results of one more integration by parts (there are no hxx terms). We 
also note that if one considered problem (P) with nonlinearity f{h) = 
\h\^ with < < 3, then Theorems [T] and [2] would hold for general 
nonnegative initial data Hq G H^{Q); no "finite entropy" assumption 
would be needed [3 [3]. Finite entropy conditions (/ hl~^ dx < oo and 
/ /iq"*"^"" dx < oo) would be needed to obtain the results for n > 3. 
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3.1 Regularized Problem 

Given S,e > 0, a regularized parabolic problem, similar to that of 
Bernis and Friedman [5], is considered: 

ht + {f5e{h){aoh^xx + aihx + a2w'{x)))^ + ash^ = 0, (3.16) 
's,e) \ g(-a,t) = g(a,t) for t > 0, 2 = M, (3.17) 
h{x,0) = ho^seix), (3.18) 



where 



fseiz) := fe{z) + 5 = ^ + 5 \/zeR\ 6>0, e>0. (3.19) 



The 5 > in (3.19) makes the problem (3.16) regular (i.e. uniformly 
parabolic). The parameter e is an approximating parameter which has 
the effect of increasing the degeneracy from f{h) ~ \hf to feih) ~ /i^. 
The nonnegative initial data, ho, is approximated via 



ho,5e = ho,5 + e C^+^(n) for some 9 G (0, 2/5) and 7 from ( |3.4| ) 

%^^(-a) = %^^(a)forz = 0;3, 
^o,(5e ho strongly in H^{il) as 5, £ — > 0. 

(3.20) 



The e term in (3.20) "lifts" the initial data so that it will be positive 
even if 5 = and the 6 is involved in smoothing the initial data from 

H\n) to c^+^(fi). 

By Eidelman [T71 Theorem 6.3, p. 302], the regularized problem has 
a unique classical solution hsg G C^'^^'^~^^^'^{^1 x [0,r5e]) for some time 
Tse > 0. For any fixed value of 6 and e, by Eidelman [17., Theorem 9.3, 
p. 316] if one can prove an uniform in time a priori bound \hss{x,t)\ < 
Ase < 00 for some longer time interval [0,T/oc,fe] (^/oc.fe > Tse) and 
for all X G f2 then Schauder-type interior estimates [13 Corollary 2, 
p. 213] imply that the solution hse can be continued in time to be in 

Although the solution hse is initially positive, there is no guarantee 
that it will remain nonnegative. The goal is to take 5 — 0, e — >^ in 
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such a way that 1) Tioc^se — > Tioc > 0, 2) the solutions hse converge 
to a (nonnegative) hmit, h, which is a generalized weak solution, and 
3) h inherits certain a priori bounds. This is done by proving various 
a priori estimates for hse that are uniform in 6 and e and hold on 
a time interval [0,T;oc] that is independent of 6 and e. As a result, 
{hse} will be a uniformly bounded and equicontinuous (in the Cl^f'^^^ 
norm) family of functions in x [0,T;oc]- Taking 6^0 will result in 
a family of functions {h^} that are classical, positive, unique solutions 
to the regularized problem with 5 = 0. Taking e — > wiU then result 
in the desired generalized weak solution h. This last step is where 
the possibility of nonunique weak solutions arise; see [3] for simple 
examples of how such constructions applied to hf = —{\h\"'hxxx)x can 
result in two different solutions arising from the same initial data. 



3.2 A priori estimates 

Our first task is to derive a priori estimates for classical solutions of 
(3.16)-(3.20). The lemmas in this section are proved in Section [A} 
We use an integral quantity based on a function Gse chosen so that 



and Gs^{z) > 0. 



(3.21) 



This is analogous to the "entropy" function first introduced by Bernis 
and Friedman [5]. 

Lemma 3.1. There exists 5q > 0, > 0, and time Tioc > such that 
if S E [0,6o), e G {0,eo), if hse is a classical solution of the problem 
(3.16)- [3. 20 ) with initial data ho^se, o-nd if h^^se satisfies (3.20) and 
is built from a nonnegative function ho that satisfies the hypotheses of 
Theorem^ then for any T E [0,Tioc] the solution hse satisfies 



I {hl^^{x,T)+^^(^^ + 26)Gse{hss{x,T))}dx (3.22) 
+ ao fse{hse)h'j^ dxdt ^ Ki < oo, 
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J Gss{hss{x, T)) dx + ao J J h]^,^^ dxdt < K2 < oo, (3.23) 



and the energy Eseit) (see (3.13)) satisfies: 



lse{T) + j j f5eihse){aohse,xxx + aihse,x + a2w'f dxdt (3.24) 



^ Co + i^sT 

where K3 = \a2a3\ \\ w'\\ooC < 00. The time Tioc and the constants Ki, 
K2, Co, and K3 are independent of 6 and e. 

The existence oi 5q, £9, T;oc, -f^i, -^^2, and is constructive; how to 
find them and what quantities determine them is shown in Section |Xj 
yields uniform-in-5-and-5 bounds for J h^^^, J Gseihse), 



Lemma 



3.1 



// ^lexx^ and J f f5e{hse)h1^^^^. However, these bounds are found in a 
different manner than in earher work for the equation ht = —{\h\"'hxxx)x, 



for example. Although the inequality (3.23) is unchanged, the inequal- 
ity (3.22) has an extra term involving Gse- In the proof, this term 



was introduced to control additional, lower-order terms. This idea 
of a "blended" ||/i:i.||2-entropy bound was first introduced by Shishkov 
and Taranets especially for long-wave stable thin film equations with 
convection pi]. 

The final a priori bound uses the following functions, parametrized 
by a, 



+ 



{a-l){a-2) ' (o-3)(q:-2)' 



(C(")(z))" 



(3.25) 



Lemma 3.2. Assume sq and Tioc are from Lemma 3.1, 6 = 0, and 



e G (0, £0)- Assume kg is a positive, classical solution of the problem 
(3.16)-(3.20) with initial data Hq^ satisfying Lemma 3.1 Fix a G 



1^—1/2, 1) with a ^ 0. If the initial data h^ ^ is built from Hq which also 
satisfies 

(3.26) 



Hq ^{x) dx < 00 
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then there exists and Ti^J with < Eq < Eq and < T/^" < Tioc 
such that 

{hl^{x,T) + Gi"\K{x,T))}dx (3.27) 
+ // [Ph'^hl^^ + -fh'^-^hl^] dxdt^Ki<oo 



Qt 

holds for all T G [0, T^^J] and some constant K4 that is determined by 
a, Eq, ao, ai, 02, w' , Q and Hq. Here, 

g^Uo if a E (0,1), ^iao^ ^/«e(0,l), 

Vo^) ^/a G (-1/2,0), ^ \ao^i±?fi^ ^/ a G (-1/2, 0). 

Furthermore, 

a+2 a+2 

he^ eL\Q,Tioc-,H\n)) and he' e L''{0,Ti,,;Wlin)) (3.28) 
with a uniform-in-E bound. 

The a-entropy, / G^f^\h) dx, was first introduced for a = —1/2 in 



[TU] and an a priori bound like that of Lemma 3.2 and regularity results 
like those of Theorem [2] were found simultaneously and independently 
in [3] and |7]. 

3.3 Proof of existence and regularity of solutions 



Bound (3.22) yields uniform L°° control for classical solutions h^e, al- 
lowing the time of existence Tioc,5e to be taken as Tioc for all 6 G (0, 5o) 
and E G {0,Eo). The existence theory starts by constructing a classi- 



cal solution hse on [0,Tioc] that satisfy the hypotheses of Lemma 3.1 
if (5 G (0, 5o) and e G (0,£o)- The regularizing parameter, 6, is taken 
to zero and one proves that there is a limit h^ and that is a gen- 
eralized weak solution. One then proves additional regularity for h^; 
specifically that it is strictly positive, classical, and unique. It then 



15 



follows that the a priori bounds given by Lemmas 3A_^ and 3^ apply 
to he- This allows us to take the approximating parameter, e, to zero 
and construct the desired generahzed weak solution of Theorems [T] and 

m 



Lemma 3.3. Assume that the initial data /lo.e satisfies (3.20) and is 
built from a nonnegative function Hq that satisfies the hypotheses of 
Theorem^ Fix 5 = and e G (0,£:o) where is from Lemma 3.1. 
Then there exists a unique, positive, classical solution on [0, T/oc] of 
problem (Po,e), see (3.16]-(3.20), with initial data /io,e where Tioc is 



the time from Lemma 3. 1 



Proof. Arguing the same way as Bernis & Friedman [5] one can con- 
struct a generalized weak solution h^. We now prove that this 
is a strictly positive, classical, uni que s olution. This uses the entropy 
/ Gse{hse) and the a priori bound (3.23). This bound is, up to the coef- 



ficient oq, identical to the a priori bound (4.17) in [5j. By construction, 
the initial data /io,£ is positive (see (3.20)), hence / (^^(/io.e) dx < oo. 
Also, by construction fs{z) ~ for z <C 1. This imphes that the gen- 
eralized weak solution kg is strictly positive [HI Theorem 4.1]. Because 
the initial data /io,e is in C^{^1), it follows that kg is a classical solution 
in C^'}{Qtioc)- The proof of Theorem 4.1 in [5] then implies that is 



unique. 



□ 



Proof of Theorem\I\ As in the proof of Lemma 3.3, following [5], there 
is a subsequence {sk} such that hsf, converges uniformly to a function 



h G C l^t'^ ^^ which is a generalized weak solution in the sense of Defi- 
nition 3.1 with f{h) = \h\'^. 

The initial data is assumed to have finite entropy: /l//io < C)0. 
This, combined with f{h) = \h\^, imphes that the generahzed weak 
solution h is nonnegative and the set of points {h = 0} in Qti^c has 
zero measure [5l Theorem 4.1]. 

To prove (|3.14l), start by taking T 



Tioc in the a priori bound 



(3.24). As Ek — > 0, the right-hand side of (3.24) is unchanged. First 
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consider the — limit of 



\oc) - 2a2w{x)hs^{x,Tioc)dx. 



the second and third terms 
0. The bound yields 



By the uniform convergence of /i^j, to h 
in the energy converge strongly as Ek — 
a uniform bound on {J^hl^ ,j.{x,Tioc) dx}. Taking a further refine- 
ment of {sk}, yields hs^,xi-,Tioc) converging weakly in L^{Q). In a 
Hilbert space, the norm of the weak limit is less than or equal to the 
liminf of the norms of the functions in the sequence, hence E,o{Ti„^) < 



loc, 

liminf(.^^o ££^(T;oc ). A u niform bou nd on// fe{K) {aoK^^^x + . . .f dx 
also follows from ( 3.24[ ). Hence \/ feAheu) {a^hst^xxx + • • • ) converges 
weakly in L'^{Qti^J, after taking a further subsequence. It suffices to 
determine the weak limit up to a set of measure zero. Because h > 
and {h = 0} has measure zero, it suffices to determine the weak limit 
on {h > 0}. 

The regularity theory for parabolic equations allows one to argue 
that h e C^'l{7), and the weak limit is /i^/^ [aohxxx + ■ ■ ■) on {h > 0}. 
Using that 1) the norm of the weak hmit is less than or equal to the 
liminf of the norms of the functions in the sequence and that 2) the 
liminf of a sum is greater than or equal to the sum of the liminf s. 



results in the desired bound (3.14). 



It follows from ( |3.23| ) that h^^^xx converges weakly to some v in 
L'^{Qtioc)} combining with strong convergence in L'^{Q,T; H^{il)) of 
hgr^, to h by Lemma B.l and with the definition of weak derivative, we 
obtain that v = hxx and h G L'^{0,Tioc; H'^i^)) that imph es (|3.12 ). 
Hence he,t ^ h weakly in L^{0,T; {H\n))') that implies i\3.6]. By 
Lemma [R2l we also have h G C{[0,Tioc\, L'^{^)). 



□ 



Proof of Theorem^ Fix a E (—1/2, 1). The initial data ho is assumed 



3.2 



holds 



to have finite entropy / G^f^\ho{x)) dx < oo, hence Lemma 
for the approximate solutions {he,,} where this sequence of approximate 
solutions is assumed to be the one at the end of the proof of Theorem [T| 
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By (pSl), 



a+2 
he,' 



and 



a+2 



h. 



is uniformly bounded in Sk in L^(0, T/oc; 



is uniformly bounded in in L'^{0,Tioc;W^{^l)). 



Taking a further subsequence in {sk}, it follows from the proof of [TSl 
Lemma 2.5, p. 330], these sequences converge weakly in L^(0, T^oc; H'^{^)) 

a+2 a+2 

and L'^{0,Tioc;Wl{Q)), to h 2 and /i 4 respectively. □ 



4 Long-time existence of solutions 

Lemma 4.1. Let h G be a nonnegative function such that 

J h{x) dx = M >0. Then 



2 4 



\n\ ■ 



Note that by taking /i to be a constant function, one finds that the 
constant M^/|fi| in (4.1) is sharp. 



Proof. Leiv = h- M/|fi|. By (A.3) 



ll,2(n) 



v1 dx 



\v\ dx 



Hence, 



1 



dx + 



1 

3 1 I Kidx ) ^(2M)i + 



n 



m ■ 



□ 
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Lemma 4.1 and the bound (3.14) are used to prove control of 



the generahzed weak solution constructed in Theorem [T] 

Lemma 4.2. Let h be the generalized solution of Theorem\I\ Then 



4 



where £o(0) is defined in (3.13), M = f Hq, K = |a2a3|||w 

\a2\\\w\\^M 



(4.2) 

-,C and 



2^6 (ao+ai)^/^ , gQ+ai 

^ 2\n\ 



if ao + 0.1 < 0, 
otherwise. 



Note that if the evolution is missing either linear or nonlinear ad- 
vection (02 = or = or 03 = 0) then Lemma 4.2 provides a 
uniform-in-time upper bound for \\h{-,Tioc)\\m. 



For the equation (1.3) which models the flow of a thin fllm of liquid 

.a 



3 ' 



a2 = 



on the outside of a rotating cylinder one has ao = ai 
03 = 1, w{x) = sinx, and = 2ti. In this case, the bound (4.2) 
becomes 



||M.,T,„,)||^.(^) ^ £o(0) + |CT,,, + fM + M^ I) 



where 2£o(0) = / (x/3 {h\^ - h^) + 2/^/3 sin(a;) ho )dx. The bound 
(|4.2[) actually holds true for aU times for which h is strictly positive. 



Recalhng the definition (1.1) of x, one sees that the control is lost 
as X — ^ (i.e. as a/^upRuj) — > 0), for example, in the zero surface 
tension limit. 



Proof By (3.13) 



f j hi{x, T) dx = £o(r) + f / h^x, T) dx + a2 / h{x, T) w{x) dx. 
a 



The hnear-in-time bound (3.14) on £o(T;oc) then imphes 



\\h{;Tioc)fHi <^o{0) + KTioc+'^ J h"" dx+MM^M. (4.3) 
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with K 



c. 



Case 1: + cii < The third term on the right-hand side of (4.3) 
is nonpositive and can be removed. The desired bound (4.2) follows 
immediately. 



Case 2: ao + cti > By Lemma 4.1 and Young's inequality 



/ dx< ( 63M3 



h^dx ] + 



\n\ 



<f y hl{x,Tioc) dx + ( 



2 f 2V6(ao+ai)3/2 qQ+qi 

2|n| 



(4.4) 



Using this in (4.3), the desired bound (4.2) follows immediately. □ 

This control in time of the generalized solution is now used to 
extend the short-time existence result of Theorem [T] to a long-time 
existence result: 



Theorem 3. Let Tg be an arbitrary positive finite number. The gen- 
eralized weak solution h of Theorem [I] can be continued in time from 
[0, Tioc] to [0, Tg] in such a way that h is also a generalized weak so- 
lution and satisfies all the bounds of Theorem [7] ( with Tioc replaced by 

Similarly, the short-time existence of strong solutions (see Theorem 
[2| can be extended to a long-time existence. 

Proof. To construct a weak solution up to time Tg, one applies the 
local existence theory iteratively, taking the solution at the final time 
of the current time interval as initial data for the next time interval. 
Introduce the times 

N-l 

= To<Ti<T2<-- - <Tn <... where := ^ Tn,ioc (4.5) 

n=0 
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and Tn^ioc is the interval of existence (A. 12) for a solution with initial 
data h{-, Tn): 



- nioc 



40C6 



minjl, i^j hl{x,T^) + 2fG^{h{x,Tn)) 



-2 ■ 



(4.6) 

The proof proceeds by contradiction. Assume there exists initial 
data /zq, satisfying the hypotheses of Theorem [TJ that results in a weak 
solution that cannot be extended arbitrarily in time: 



Tn,loc — T* < 



oo 



lim Tnjoc = 0. 



k=0 



From the definition (4.6) of T„_;oc, this implies 



lim [ihl{x,T^) + 2fGoih{x,T„)))dx = cx). (4.7) 
By (|42( and ^J§, f j hl{x,Tn) dx < E^{Tn-i)+K Tn-i,ioc+K^. 



£o(^n-l) ^ £o(^n-2) + K Tn-2,loc- 



Combining these, 



f j hl{x, Tn) dx < £o(T„_2) + K {Tn-2,loc + T, 



n—l.loc. 



Continuing in this way, f j hl{x, Tn) dx < £o(0) + KTn + K3. 

(4.8) 

By assumption, T„ — >^ T* < 00 as n — 00 hence / hl{x, Tn) dx remains 
bounded. Assumption (4.7) then implies that / Go{h{x,Tn)) dx 00 
as n — > 00. 
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To continue, return to the approximate solutions h^. By (A.8) 

j GsiK{x,Tn,e)) dx < j Ge{hsix,Tn-i,s)) dx 



(4.9) 



+ C5 J max<l,y" T) c/x > dT 



Tn-l,e 



Using (3.24), one proves the analogue of (4.2) for all T G [0,T/oc,e] and 
then the analogue of (4.8) for all T G [0,T„^e]. Using this bound, 

j j hl^i^, T) dxdT < ^ j £,(0) + KT + K^dT 



4_ 



[£,(0) + + f {Tn-l,e + r„,,)] Tn-l,loc,e. (4.10) 



Replacing K^^ by a larger value if necessary and using (4.10) in (4.9), 

Ge{hs{x,Tn,e)) dx 



< J GeiK{x, Tn-i^s)) dx+{a + P{Tn-i,e + Tn^e)) ^n- 



(4.11) 



lloc.e 



for some a and [3 which are fixed values that depend on the co- 
efficients of the PDE, and (possibly) on the initial data /lo.e- Taking 
efc — s> in the sequence {sk} that was used to construct h yields 

j Go{h{x,Tn))dx < J Goih{x,Tn-l))dx+{a + P(Tn-l +Tn)) Tn-l,loc- 

n a 

(4.12) 

Applying (4.12) iteratively and using that <T*^ 

j Go{h{x,Tn)) dx < J Go{ho{x)) dx+ {a + P2T*)Tn. (4.13) 
n n 

Hence / Go{h{x,Tn))dx < oo as n — »• cxo, finishing the proof. □ 
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Under certain conditions, a bound closely related to (4.2) implies 
that if the solution of Theorem [T] is initially constant then it will remain 
constant for all time: 



Theorem 4. Assume the coefficients ai and 02 in (1.8) satisfy ai > 0, 
a2 = and \^\ < 4ao/|ai|. // the initial data is constant, Hq = C > 0, 
then the solution of Theorem^ satisfi,es h{x,t) = C for all x E Cl and 
all t > 0. 

The hypotheses of Theorem |4] correspond to the equation is long- 
wave unstable (ai > 0), there is no nonlinear advection (02 = 0), and 
the domain is not "too large". 

Proof. Consider the approximate solution he. The definition of £e(T) 
combined with the linear-in-time bound (3.24) implies 



J hl^{x,T) dx < Ee{0) + KT+^-^ j hldx+\a2\\\w\\^M, (4.14) 



where = Jho^^dx. Applying Poincare's inequality (A. 2) to 

hs - Mj\n\ and using J h^dx = J v"^ dx + M^/\n\ yields 



2 8 



j h%{x,t) dx < ^e{0) + KT,,i,,+ \-^ + \a2\\\w\\^M,. 



If /io,e = Ce = C + and 02 = (hence K = 0) this becomes 

jhUx^T) dx<ia,-\a,\)^. 



ao 
2 



If ai > and < 4ao/ai then / hl^{x, T) dx = for all T e [0, T^^ioc] 
and that this, combined with the continuity in space and time of h^, 
imphes that /i^ = on Qt^j^^. Taking the sequence {e^} that yields 
convergence to the solution h of Theorem [Tj h = C on Qti^^. 

□ 
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5 Strong positivity of solutions 



Proposition 5.1. Assume the initial data Hq satisfies Hq^x) > for 
all X & u C Q where u is an open interval. Then the weak solution h 
from Theorem [I] satisfies: 

1) h{x,T) > for almost every x E uj, for all T G [0,T;oc]; 

2) h{x, T) > for all x E uj, for almost every T G [0, TjoJ- 



The proof of Proposition |5.1| depends on a local version of the a 



priori bound (3.23) of Lemma 3.1 



Lemma 5.1. Let u (1 Q be an open interval and ( G C^(fi) such that 
( > on u, supp( = uj, and (C^)' = on dQ. If u = fl, choose ( such 
that C{-a) = C(a) > 0. Let ^ := C^. 

// the initial data Hq and the time Tioc are as in Theorem [I] then for 
all T G [0,T;oc] the weak solution h from Theorem^ satisfies 



e(x) 



1 



h{x,T) 



dx < OO 



(5.1) 



n 



The proof of Lemma 5T is given in Appendix|X] The proof of Propo- 
sition [5T] is essentially a combination of the proofs of Corollary 4.5 and 
Theorem 6.1 in |5] and is provided here for the reader's convenience. 



Proof of Proposition \5.1\ Choose the locahzing function ((x) to satisfy 



the hypotheses of Lemma 5.1 Hence, (5.1 ) holds for every T G [0, Tioc]- 
First, we prove h{x,T) > for almost every x G cu, for all T G 
[0,T;oc]. Assume not. Then there is a time T G [0,T;oc] such that the 
set {x I h{x,T) = 0} rio; has positive measure. Then 



oo > 



e(x) 



h{x,T) 



dx > 



^x) 



1 



h{x,T) 



dx 



oo. 



n 



{h{-,T)=o}nuj 



This contradiction implies there can be no time at which h vanishes 
on a set of positive measure in u, as desired. 

Now, we prove h{x,T) > for all x G cu, for almost every T G 



[0,Tioc]. By (3.12), h^^{-,T) G L^{n) for almost all T G [0,Tioc] hence 
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h{-,T) E C3/2(fi) for almost all T G [0,Tioc]. Assume Tq is such that 
/?.(■, To) e C^/^(fi) and /i(xo,To) = at some xq G tu. Then there is a 
L such that 

h{x,To) = \h{x,TQ) - h{xo,To)\ < L\x - xq\^^'^. 

Hence 

This contradiction imphes there can be no point xq such that h{xQ, Tq) = 
0, as desired. Note that we used ^ > on a; and G a; to conclude 
that the integral diverges. 

□ 

We close our discussion by illustrations of positivity and long time 
existence via numerical simulations of the initial value problem for 
different regimes of the PDE. 

Figure [2] considers the PDE with no advection, ht + [h^ih^xx + 
^&hx))x = 0. The PDE is translation invariant in x and constant 
steady states are linearly unstable. As a result, any non-constant be- 
haviour observed in a solution starting from constant initial data would 
be due to growth of round-off error. For this reason, non-constant ini- 
tial data is chosen: ho{x) = 0.3 + 0.02 cos(x) + 0.02 cos(2x). The 
and norms of the resulting solution appear to be converging to hm- 
iting values as time passes and long-time hmit of the solution appears 
to be four steady-state droplets of the form acos(4x + (j)) + b for ap- 
propriate values of a, 0, and b. Like the PDE, the simulation shown 
respects the symmetry about x = of the initial data. 

Figure [3] shows the evolution from constant initial data for the PDE 
with nonlinear advection but no linear advection: ht+{h^{hxxx+^G hx — 
8 cos(x)))a; = 0. The long-time hmit appears to be a steady state which 
is zero (or nearly zero on [— vr, 0] ) with a droplet supported within (0, vr) 
and centred roughly about the mid point (x = it/2). 
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Figure 2: The evolution equation with no linear or nonlinear advection, 
ht + {h^ihxxx + 16 /ii:))^ = 0, corresponding to Oq = 1, ai = 16, and 
«2 = 03 = 0. The initial data is /lo (3;) = 0.3+0.02 cos(x)+0.02 cos(2x). 
Left plot: the solution at times t = (dashed line), t = 12, 12.5, 13, 15 
(solid lines), and t = 140 (heavy line). Right plot: the Li^ and 
norms plotted function of time. 



Finally, Figure |4] shows the evolution resulting from the same con- 
stant initial data for the PDE with both linear and nonlinear advec- 
tion: ht + {h^ihxxx + 16 /ix — 8 cos{x)))x + ^hx = 0. The long-time limit 
appears to be a strictly positive steady state. 

We close by noting that the PDE considered in Figure [4] corresponds 
to coefficient 03 = 3 in the PDE (|1.8l). As we increase the value of 03 



we find there appears to be a critical value past which the solution 
appears to converge to a time-periodic behaviour rather than a steady 
state. 



A Proofs of A Priori Estimates 



The first observation is that the periodic boundary conditions imply 



that classical solutions of equation (3.16) conserve mass: 



hse{x, t) dx 



hf) seix) = ^fiF < 00 for all t > 0. 



(A.l) 



n 
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Figure 3: The evolution equation with nonhnear advection but no 
linear advection, ht + {h?{hxxx + 16 /i^ — 8 coii^x)))^ = 0, corresponding 
to 0,0 = 1, fli = 16, 02 = 8, and 03 = 0. The initial data is ho^x) = 0.3. 
Left plot: the solution at times t = (dashed line), t = 0.5,1,2,10 
(solid lines), and t = 3000 (heavy line). Right plot: the and 
norms plotted function of time. 



2 



1.5 



0.5 








Figure 4: The evolution equation with botlTlinear and nonlinear ad- 
vection, ht + {h?{hxxx + 16 /la; — 8cos(a;)))a; + ?>hx = 0, corresponding 
to oo = 1, oi = 16, 02 = 8, and 03 = 3. The initial data is ho{x) = 0.3. 
Left plot: the solution at times t — (dashed line), t = 0.5,1,2,4 
(solid lines), and t = 20 (heavy line). Right plot: the and 
norms plotted function of time. 
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Further, (3.20) implies M^g ^ M = J ho as e,5 ^ 0. The initial data 
in this article have M > 0, hence Mse > for 5 and e sufficiently small. 

Also, we will relate the norm of h to the norm of its zero-mean 
part as follows: 



\h{x)\< h{x)-^^ 



+ 



M_ 



p-i 



MP 



where v := h — M/\il\ and we have assumed that M > 0. We will 
use the Poincare inequality which holds for any zero-mean function in 

||^||^<&i||t':.||^ l<p<cx) (A.2) 

with bi = \n\P/{p2P-'^). 

Also used will be an interpolation inequality [22l Th. 2.2, p. 62] for 
functions of zero mean in H^{Q): 



\v\\P < bo \\v. 



I ip 



\(l-a)p 



(A.3) 



where r > 1, p > r, 

« = ^, b2 = {l + r/2y 
It follows that for any zero-mean function v in H^{fl) 



\h\\l<b4K\\l + b5Ml 



5e 



(A.4) 




6i |fi|(2-p)/p if l<p<2 

^(p+2)/2 if 2 < p < CX) 



2P-' 6, 



1^1 



p-1 



To see that (A.4) holds, consider two cases. If 1 < p < 2, then by 
(A.2), ||t;||p is controlled by ||fxilp- By the Hold er inequality, || 

^x||p is 

then controlled by Hf^lh- If p > 2 then by (A.3), ||f ||p is controlled by 



|l-a 



■'x\\2\ 



where a = 1/2 — 1/p. By the Poincare inequality, ||w 



is controlled by ||w 



|l-a 



x\\2 
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Proof of Lemma \3.1\ In the following, we denote the classical solution 
/ife by h whenever there is no chance of confusion. 



To prove the bound (3.22) one starts by multiplying (3.16) by —h 



integrating over Qt, and using the periodic boundary conditions (3.17) 
yields 

I J hl{x,T) dx + ao JJ fs£{h)hl^^dxdt (A. 5) 

n Qt 

= I y ho^se,x^ix) dx - ai j j feih)hxha:xx dxdt + Sai j j hl^ dxdt 

n Qt Qt 

- a. jj Mhyh... dxdt - Sa, If u/ ft„ d.dt. 

Qt Qt 



By Cauchy and Young inequalities, due to (A.2)-(A.4), it follows 



from (A. 5) that 



1 1 hUx,T)dx + f II f 



'&e{h)hl^^ dxdt 



— \ I ^o,fe,a;^ dx + C3 II hl^ dxdt + C4 y max < 1, | y dx 



(A.6) 



dt 



where ci = bj/S + 64/2, C2 = Ml h/2, C3 = + 5|ai| 
C4 = ^^ci + ^Ww'Wlb^ + + ^\\w'\\lb,Ml + 6^\\w'\ 

^ ao ao " , "1 " "0 " 1100 o OS ' ao II I 

Now, multiplying (3.16) by G'^^{h), integrating over Q^, and using 



12 

l2- 



the periodic boundary conditions (3.17), we obtain 



I Gseih(x,T)) dx+ao H h^^dxdt = I Gse{hQ,5e)dx+ai H h^^dxdt 

n Qt ^ Qt 

— a-i {Gseih))^ dxdt + a2 / / w'h^dxdt. (A. 7) 
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By the periodic boundary conditions, we deduce 

Gseih{x, T)) dx + ao J J h^^ dxdt 
n Qt 



^ J Gseiho^se) dx + C5 J max < 1, J hl{x,t) dx\ dt, {A.i 
n K n ) 



where C5 = |ai| + |a2| ||w^'||2- Further, from (A. 6) and (A.8) we find 



J hldx+^ j G&e{h) dx + ao jj fse{h)hl^^dxdt < j ho^se,x^dx 



n 

3- 



+ ^ J G5eiho,se)dx + ^6 j max S 1' ( y hl{x,t) dx\ \ dt (A. 9) 

Q I \q / j 

where cq = 20305/00 + 2c4. Applying the nonlinear Gronwall lemma 
[12] to 

T 

v{T) < v{0) + ce / max{l,v%t)} dt 







with v{t) = J{hl{x,t) + 2c3/ao Gse{h{x,t))) dx yields 

j hl{x,t) +2f^G&e{h{x,t)) dx (A.IO) 

< v^max |l,y"(Ve,x^(a;) + 2^G5,(VeW)) rfxj = fsT^^ < 00 



for all t e [0, T^e^/oc] where 



Tse,ioc '■= 4^ min | ^/ ^^o,Se/{x) + 2^Gse{ho^5e{x))) dx^ I 



(A.11; 
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Using the S ^ 0,e ^ convergence of the initial data and the choice of 



6 G (0,2/5) (see (3.20)) as well as the assumption that the initial data 
ho has finite entropy (3.11 ), the times Tse,ioc converge to a positive limit 



and the upper bound K in (A. 10) can be taken finite and independent 



of S and e for 6 and e sufficiently small. (We refer the reader to the end 



of the proof of Lemma 5.1 in this Appendix for a fuller explanation of 



a similar case.) Therefore there exists 5o > and Eq > and K such 



that the bound (A. 10) holds for all < S < Sq and < e < Eq with K 



replacing Kse and for all 



0<t<Tia 



9_ 
10 



lim Tse,ic 



(A.12) 



Using the uniform bound on J hi that (A. 10) provides, one can find 
a uniform-in-5-and-e bound for the right-hand-side of (A. 9) yielding 



the desired a priori bound (3.22). Similarly, one can find a uniform- 



in-5-and-e bound for the right-hand-side of (A.8) yielding the desired 



a priori bound (3.23). 



To prove the bound (3.24), multiply (3.16) by —a^h^x — aih — 



integrate over Qt, integrate by parts, use the periodic boundary con- 



ditions (3.17), and use the mass conservation (see (A.l)) to find 



£fc(T) + jj fse{h){aohxxx + ^ihx + a2w'{x)f dxdt 

Qt 

< £fe(0) + laaaal Ik'lloo + 2m) T. (A.13) 



Hence the desired bound (3.24) is obtained if the constant 

K3= \a2a3\\\w'\\UM^^/K~i + 2M). 

The time Tioc and the constants Ki, K2, and are determined by Sq, 

Eq, Oo, di, 02, w', and ho. □ 



Proof of Lemma \3^ In the following, we denote the positive, classical 
solution hs by h whenever there is no chance of confusion. 
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Multiplying (|3.16[) by {Gi°'\h)y, integrating over Qt, taking 6^0, 



and using the periodic boundary conditions (3.17), yields 

Gi''\h{x, T)) dx + aoJJ h'^hl^ dxdt + 



Qt 



Qt 



J Gi^'Xho,) dx + ai Jj h^hl dxdt jj h'^+^w" dxdt. 



(A.14) 



Qt 



Case 1: < a < 1. The coefficient multiplying ff in (A.14) 

is positive and can therefore be used to control the term J f h h'^ on 
the right-hand side of (|A.14l). Specifically, using the Cauchy- Schwartz 



inequality and the Cauchy inequality, 

aoa(l-a) / / ,a-2v4 



ai / / h'^hi dxdt ^ ^^2i!ip^Z / / h^--^hl dxdt+ 1 1 dxdt. 



(A.15) 



Using the bound (A.15) in (A.14) yields 



j Gi"\hix,T))dx + ao jj h"hl^dxdt + ao^^^ jj h'^-^hldxdt 



Qt 



(A.16) 



< j Gi"\hos) dx + ^^^^^ j I' h-^^ dxdt + »^ j j h-^'dxdt. 



By (A.4) 



Gi^'Xhix, T)) dx + ao / / h^'hl^ dxdt + ao^^^ / / h^'-^hl dxdt 



WT 



< j G'f\hoe) dx + dij max i 1, j ^ /i^ c/x j \ dt (A. 17) 



32 



where 



di = 64 



3a? L |a2|||«."||. 



+ 



2aoa(l— «) 1+a 



3a? 



2aoo(l— a) 



+ 2 _|_ |a2|H«'"|joo + l 

£ 1+a £ 



Using the Cauchy inequahty in (A. 9) and taking 5 — yields 



h^. dx + Oo J J fe{h)hl^^ dxdt 



(A.18) 



Qt 

< I hi dx + 



2a? 
ao 



a; ao 



Qt 



Qt 



Applying the Cauchy-Schwartz inequality and (A.4) yields 



hi dx + ao J J fe{h)hl^^ dxdt < J h^^,^ dx 



Qt 



+ 



aoa(l— q) 



j j h^ ^/i^ dxdt + d2 J max llAJ h^ dx 



dt 



Qi 



where 



6af 2a| n / ii2 



d2 = bd { -3-Tr — V + — \\w' 



\ aga(l-a) ^'^e 



Adding / Gi''\h{x,T)) to both si des of ( |A.18| , ao JJ h'^hl^ to the re- 
sulting righthand side, and using (A. 17), 



y /i^(x,T)dx + y G^^\h{x,T))dx + ao jj fe{h)hl,^dxdt (A.19) 



Qt 



T 



— J f^oe,xdx + J G^^\hQ^) dx + ds J max < 1, I y hi dx 
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where = di + d2- Applying the nonlinear Gronwall lemma [H] to 



w(T) < v{0) + d3 J max{l, v^-^'^^t)} dt 



with v{T) = j\hl{x,T) + Gi"\h{x,T))) dx yields 



{hl{x,T)+G^^\h{x,mdx 



(A.20) 



< 46-- max 1 1, y {ho,sl{x) + G(")( c^^j = < oo 



for all T: 



n < T < T*-"-* — 



^minjl, f / {ho,el{x)+Gi"\ho,s{x))) dx 



6— a 
2 



The bound (A.20) holds for all < e < £o where Eq is from Lemma 3.1 



and for all t < min{T;oc, T^°ilc} where Tioc is from Lemma 

Using the e convergence of the initial data and the choice of 6* G 



3.1 



(0,2/5) (see (3.20)) as we ll as the assumption that the initial data ho 
has finite a-entropy (3.26), the times T^'^l^ converge to a positive limit 
and the upper bound in (A.20) can be taken finite and independent 



of e. (We refer the reader to the end of the proof of Lemma 5.1 in this 
Appendix for a fuller explanation of a simi. 



exists eo"-* and K such that the bound (|A.20|) holds for all < e < e, 
with K replacing and for all 



ar case. , 



Therefore there 

-0 



<t<T, 



(a) 



loc 



mm 



T, 



loc 



' 10 



limT 



(a) 



e.loc 



(A.21) 



where Tioc is the time from Lemma 3.1 Also, without loss of generality. 



^o"^ can be taken to be less than or equal to the Eq from Lemma 



3.1 
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Using the uniform bound on f h^. that (A. 20) provides, one can find 



a uniform-in-e bound for the right-hand-side of (A. 17) yielding the 
desired bound 



G^^\h{x,T))dx+aQ il K^hl^dxdt+aQ 



a{l—a) 
6 



h^'-'^ht dxdt < Ki 



Qt 



Qt 



which holds for all < e < ej, "^ and all < T < T; 



(A.22) 



It remains to argue that (A.22) imphes that for all < £ < e, 



(a) 



that he^'^'^^ and h^^^'^^^'^ are contained in balls in L^(0, T; if^(fi)) and 
L^(0,T; Wl{il)) respectively. It suffices to show that 



Qt 



dxdt < K, 



^^a/4+l/2^)4 ^^^^ ^ ^ 



for some K that is independent of £ and T. The integral // (/ie^^^^)^ 
is a hnear combination of// // h^'^hlh^x^ and // h°'hl^. In- 

tegration by parts and the periodic boundary conditions imply 



l-a 
3 



h^'-'^ht dxdt 



h 



Qi 



Hence //(/i" is a linear combination of // /i" ^/i^, and fjh'^hl^. 



(A.23) 



a/2+1.2 
Jxx 

By (A.22), the two integrals are uniformly bounded independ ent o f e 
and T hence // {he^'^^^)l^ is as well, yielding the first part of (3.28). 



The uniform bound of JJ{he^^~^^^'^)'^ follows immediately from the 



uniform bound of JJ ^/i^, yielding the second part of (3.28). 



Case 2: — ^ < a < 0. For a < the coefficient multiply- 



ing ///i" '^h^ in (A. 14) is negative. However, we will show that if 



a > —1/2 then one can replace this coefficient with a positive coeffi- 
cient while also controlling the term / j h'^h'^ on the right-hand side of 
dATl). 



Apply ing the Cauchy-Schwartz inequality to the right-hand side 
of (A.23), dividing by J JJ h°'~'^h'^, and squaring both sides of the 



35 



resulting inequality yields 



h'^-'^ht dxdt < 



(l-a)2 



h'^hi^dxdt Va<l. (A. 24) 



Using (A.24) in (A.14) yields 



Gi''\h{x,T))dx + ao^ jjh^hl^dxdt (A.25) 

Qt 

< j G^^\ho,) dx + ai j j h^hldxdt+\^^\\w"\\^ ff K'+^dxdt. 



Note that if a > —1/2 then all the terms on the left-hand side of (A.25) 
are positive. We now control the term JJ h"h\ on the right-hand side 
of (|X25l). 



By integration by parts and the periodic boundary conditions 

h'^hl dxdt = -j^ h^'+^K^ dxdt. (A. 26) 



Applying the Cauchy inequality to (A. 26) yields 



ai y / h^h^ dxdt < 

Qt 



a()(l+2a) , g, 2 , 
l-a) "' "'XX ' 



2(l-a) 



2ao{l+2a)(l+a)2 



Using inequality (A. 27) in (A.25) yields 



G'f\h{x,T))dx + ao^ff^) II h^hl^dxdt 



) dxdt. 
(A.27) 

(A.28) 



< / G^:'\hoe)dx + 



AQ+2 I 1^21 va+l 



2a„(l+2a)(l+a)' 



Qt 



Adding 



ao(l+2a)(l-a) / / ^.-2^4 ^^^^ 



36 



to both sides of (|A.28p and using the inequahty (|A.24|) yields 

(1+2q!) 
'4(1-0) 

n Qt 

ao(l+2a)(l-a) / / r^a-2ii 



Gi"\h{x,T))dx + ao^j^ II h'^htdxdt 



(A.29) 



+ 



+ 



36 



j I h^'-^hl dxdt < j G^^\hoe) dx 



Qt 



2ao(l+2a)(l+a)2 



a+l I 



Qt 



Using (A.29) and (A. 4) yields 



G^^\h{x,T))dx + 



-J^d^h^hl, + '^0(l+2»)(l-a) ^.-2^4 I ^^^^ 



36 



Qt 



r+1 



< j G'f\hoe) dx + ei j max J 1, j y /i^ c/x j \ dt (A.30) 



where d = 64 (s^^^^feySW + SlU^"!!' 
+ ^l|w"||oo iVf^f+i). Recall the bound \A.18\ . As before, by the 
Cauchy inequality, 

i /" /" h^'hl dxdt < °0(l+2«)(l-») /" f ^a-2^4 ^^^^ (^ 3^^ 



&5( 



0+2 



2ao(l+2a){l+a)2 ^"e 



2af 
ao 



Qt 



Qt 



+ 



36af 



a^(l+2a)(l-a) 



Using (A.31) in (A.18) yields 



hldx + ao JJ fs{h)hl^^dxdt < J hl^ ,^dx 



Qt 



+ °o(i+2a)(i-a) 1 1 ^«-2^4 ^^^^ ^ ^ maxjl, ( I hi dx 



dt 



Qt 
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where 62 = &4 



1 \ ^\\w ''^ 



ag(l+2a)(l-a) ' ao lloo 



2a| 
ao 



\W 



'11^ M3). Just as (A.17) and (A.18) yielded (A.19), (A.30) com 



bined with the above inequality yielas 



hl{x,T)dx + j Gi"\h{x,T))dx + ao jj fe{h)hl^^dxdt (A.32) 



Q2 



< / hl^ .^dx+ / ^^"^(/loe) (ix + 63 / max< 1, ( / hi dx 



n 



where 63 = 61 + 62- The rest of the proof now continues as in the 
< a < 1 case. Specifically, one finds a bound 



{hl{x,T) + G^^\h{x,T))) dx 



(A.33) 



< 46-- max |l, j {ho,e/{x) + G(")(/^o,.(x))) cia;| = iT, < 00 



for all T: 



< T < T 



e,loc ■ 63(6-0) 



^^^minjl, ( / {ho,e/{x)+Gi"\hoAx))) dx 



The time t/^"^ is defined as in (A. 21) and the uniform bound (A.33) 



used to bound the right hand side of (A.30) yields the desired bound 



J Gi'^\h{x,T))dx+^-^^ IJh-hldxdt 



+ 



ao(l+2a)(l-a) / / 4 



36 



h'^-'hldxdt < K2. 



(A.34) 



□ 
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Proof of Lemma\5l\ In the following, we denote the positive, classical 



solution he constructed in Lemma 3.3 by h (whenever there is no chance 
of confusion). 



Recall the entropy function G'fe(z) defined by (3.21). Multiplying 



(3.16) by ^{x)G'g^{hs£), taking 5^0, and integrating over Qt yields 



^{x)Ge{h{x,T))dx — J C,{x)Ge{ho^e)dx = —as J J ^{x)G'^{h)hxdxdt 
n n Qt 

+ jj fe{h){aohxxx + aih^ + a2w'){^'G'^{h) + ^G'^{h)hx)dxdt 

Qt 

= ^3 / / i'Gs{h)dxdt + II ^'fe{h)G'^{h){aoKxx + 0'ihx + a2w')dxdt 



+ ^KiaoKxx + aiK + a2w') dxdt =: h + I2 + h- (A.35) 



One easily finds that for all £ > and all z > 

IfMG'Ml^lz, \f:{z)G'M\^2, 



j fe{s)G',{s) ds <\z^ + liiQ<e< (v^ - 3)/4 



Using these bounds, and recalhng ^ = C^, we bound I/2I: 



\h\ <jj if Chi + 71 [c'c' + c'lCx..! + c + eCx] {h' + hi)) dxdt 

Qt 

+ 2\a2\\\w'\\oo Jj C^\Cx\h dxdt + l\ai\ Jj \C\ dxdt (A.36) 

Qt Qt 

where 71 = max{102ao, 6|ai|} and < £ < {y/33 — 3)/4. Now, inte- 
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grating by parts in 1^, we deduce 



h + ao ihl, dxdt < 72 / / [C'C + C'lCxxl + C] hi dxdt 



+ 4|a2|(lk'l|oo + lk"l 



If {em+C')hdxdt (A.37) 



where 72 = max{6ao, |ai|}- Using bounds (A. 36) and (A.37) we obtain 
that 

J iG,{h,{x,T))dx ^ j ^G,ihoe)dx + C (A.38) 



where C > is independent of e > 0. Using the fact that 9 was chosen 
so that 9 <2/5 < 1/2, we have \^{x) Ge{ho,{x))\ < ^{x){Go{ho{x)) + 
c) < G{Go{hQ{x)) + c) almost everywhere in x and for all e < Eq. 
To finish the proof we apply Fatou's lemma to the left-hand side and 
Lebesgue Dominated Convergence Theorem to the right-hand side of 
(|X38|. □ 



B Results used from functional analysis 

Lemma B.l. ([2^J) Suppose that X, Y, and Z are Banach spaces, 
X (^Y C Z, and X and Z are reflexive. Then the embedding {u G 
LPo(0,T; X) : dtU G Lp^{0,T; Z),l < pi < 00, i = 0,1} m LP"(0,T;y) 
is compact. 

Lemma B.2. ( ^32[j Suppose that X, Y, and Z are Banach spaces and 
X<mYcZ. Then the embedding {u G L°°(0,T;X) : dtU G Lp{0,T;Z), 
p > 1} (E C(0, T; y) is compact. 
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